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Abstract: In this paper an initial value investigation has been made for a two dimensional MHD flow of a 
second order Rivlin-Ericksen type visco-elastic fluid through a porous medium in a parallel plate channel under 
externally applied boundary acceleration. Ihe exact solutions for the velocity, accelerations and shear stress 
have been obtained using separations of variables method. Their behaviour for different governing parameters 
has been discussed computationally based on the available physiological data. 
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I. INTRODUCTION 

External acceleration of large amplitudes imparted to the human body cause serious problems in the 
cardiovascular system, causing the impairment of certain physiological functions. Head ache, increase of 
pulse rate, loss of vision, venous pooling of blood in the lower extremities, congestion of liver and lungs 
hemorrhages in the face and neck, eyes, lungs and brain are some of the complications as a result of 
external accelerations [1, 3 & 4]. On the other hand, sudden accelerations are deliberately imparted to the body 
to obtain some beneficial effects. For example, body accelerations synchronous with heart beat, a proposed 
new method for assisting the failing heart, has been employed by Verdouw et.al [8] to reactive the failing 
heart. Experimental results obtained by Arntzentus et.al [2] on heart action in pigs indicate that blood pressure 
and cardiac output are raised when body accelerations synchronous with heart beat was applied in the forward 
direction in early systole. One of the channel walls in subjected to externally applied boundary acceleration. 
Against from this, blood is in constant circulation with in the blood vessels throughout the body. This 
circulatory system is the transport system of the body consisting of four type's viz., systematic circulation, the 
pulmonary circulation, the portal circulation and the coronary circulation. In these circulatory systems, the blood 
ultimately flows through a network of capillaries with porous wall having distinct functions. These capillaries 
act as loading or unloading stations in a transportation system with main function of nourishing the tissues in 
general. The flow through such net work of capillaries may be idealized as a two dimensional flow through a 
porous medium. The interaction of body accelerations with blood flow is desirable for controlling the ill effects 
of such body acceleration over blood circulation Lachlan [6]. Krishna and Rao [5] studied an initial value 
investigation has been made for a two dimensional flow of a second order Rivlin-Ericksen type visco-elastic 
fluid and the exact solutions for the velocity, acceleration, shear stress and volume flow rate have been obtained 
using transform method. Numerical computations have been carried out to discuss their behaviour for different 
parameters based on the available physiological data. Sud. V.K., et al [7] analyzed the blood flow in large and 
small arteries under the influence of externally applied periodic oscillations. The solutions for the blood 
velocity, fluid acceleration and shear stress are obtained. It has been show that high blood velocity and large 
shear stress are produced in large arteries. However, in the case of small arteries the flow is not disturbed. 
Recently Veera Krishna et.al [9] studied an initial value investigation and has been made for a two dimensional 
flow of a second order Rivlin-Ericksen type visco-elastic fluid in a parallel plate channel under externally 
applied boundary acceleration. Suneetha et.al [10] discussed the steady hydro magnetic flow of a couple stress 
fluid in a parallel plate channel bounded on one side by a porous bed under the influence of a transverse 
magnetic field and periodic body acceleration. In this paper an initial value investigation has been made for a 
two dimensional MHD flow of a second order Rivlin-Ericksen type visco-elastic fluid through a porous medium 
in a parallel plate channel under externally applied boundary acceleration. 

II. Formulation and Solution of the Problem 

We consider an incompressible viscous and electrically conducting two dimensional flow of a second 
order Rivlin-Ericksen type visco-elastic fluid in a parallel plate channel bounded by a loosely packed porous 
medium. The fluid is driven by a uniform pressure gradient parallel to the channel plates and the entire flow 
field is subjected to a uniform inclined magnetic field. Choose a Cartesian system (?, y) with boundary walls 
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y = 0 and h. The flow is uni-directional in view of the fact the pressure gradient and body acceleration are in the 
same direction. The governing equation of the motion in the non-dimensional flow with respect to the frame of 
reference is, 



du 1 dp 8 2 u a. 8 3 u ou 2 H 2 v 

— = - + v — - + — — — — u u 

dt p dx dy' p dy dt p k 

With boundary conditions 



g.e 



at 



y = h 



(2.1) 



(2.2) 



dt 

u=0 at y = 0 (2.3) 

Equation (2.1) is to be solved subjected to the conditions (2.2) and (2.3), In view of Equation (2.2), we take 

Jnt 



u=v(y).e" 
And the applied pressure gradient 



= A.e" 



dp 
dx 

Making use of (2.4) and (2.5) we solve (2.1) using the boundary conditions, 

g 



(2.4) 
(2.5) 

(2.6) 



v = — Sinnt on y = h 
n 

v = 0 on }> = 0 (2.7) 

Where, u is the velocity of the fluid, p is the density of fluid, /u is the coefficient of viscosity of the 

fluid, p is the pressure, A is an amplitude of pressure gradient, a, the coefficient of visco-elasticity g is the 

body acceleration (or) amplitude of the applied acceleration and V its angular frequency in rad. sec" 1 
From (2.1) and (2.4) we obtain 



d 2 v(y) 



M, 



-v(y)= 



dy 1 ju + a^in 
Where, = pin + ofi 2 e H 2 +- 



M 

k 

The general solution of Equation (2.8) is 



(2.8) 



v(y) = C, cosh- 



-y + C 2 sinh 



yM, __a_ 

^fju + ajn Af, 



y] ju + a^n 

Now we find the arbitrary constants C l and C 2 by using equations (2.6) and (2.7) we get 



(2.9) 



v(y) = cosh— = 



M, 



1 - cosh 



y + 



h 



jU + a { in 



H — sinn? 
n 



sinh 



M, 



. _jMy A 

. sinh , y 

yjju + ayin My 



in 



Equation (2.10) is substituted in equation (2.4) we get 



(2.10) 



— cosh — j=*= y + - 
Mj yjju + a x in 



A 
Mi 



{ 



1 - cosh 



M, 



V 



ju + ajn 



■ing 



sinh 



\M X 



.sinh 



■J/u + aJn 



A 



i 



ju + a x in 



Taking the real part of the equation (2.1 1) obtained as 



U = Aj + W, + U 2 + 



(2.11) 
(2.12) 



where 

Aj = - M 2 s'mnt 
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Aj = sinh 2 aA.cos 2 b/z + cosh 2 a/z.sin 2 b/z 
u { = M 2 sinhay. sinby. comt 
u 2 = M 2 cos]ray.cosby.sirm? 

u } = M 2 [ sinha/z.cosb/z. cosha^.sinbj-cosha/z.sinby.sinhay.cosby - 

cosha/z . co s 2 b/z . sinha/z . co shay . sinb y+ 

+cosh 2 a/z.cosb/z.sinb/z.sinhay.cosb); - sinh 2 a/z.sinb);.cosb/z.sinhaj.cosbj- 
-cosha/z.sinha/z.sin 2 b/z.cosha);.sinb);.]cosn? - ng.sinha/z.cosb/z. 
coslray. sinby. cosn? + ng.coshah.smbh.s'mhay.cosby.cosnt 

u 4 = M 2 [sinha/z.cosb/z.sinha);.cosby+ 

+cosha/z.sinb/z.sinha);.sinb);-cosha/z.cos 2 b/z.sinha/z.sinha);.cosb);- 
- cosh 2 a/z.sinb/z.cosb/z.cosha);.sinb);+sinh 2 a/z.sinb/z.cosb/z.cosha);.sinb);- 
-coshah.smhah.sm 2 bh.smhay.cosby]s'mnt - ng.sinha/z.cosb/z.sinha); 
cosbj.sinnr-ng.cosha/z.sinb/z.coshaj.sinbj.sinn? 

k 



M 2 = 



2 2 , 

p n + 



k j 



a = 



b = 



J2 



4~2 



R(S - (M 2 + D ')) + VR'CS-tlVf+D-'Xt'+l 



VR 2 (S-(M 2 +D')) 2 +lf 
VR 2 (S-(M 2 +D 4 )) 2 +1 - R(S - (M 2 + D ')) 



VR 2 (S-(M 2 +D')) 2 +lf 



Here, the non-dimensional parameters R, S, p , M, D" 1 and G are given by 



R = 



h 2 np 



ft 



is the Reynolds number, 



S = — — is the visco-elastic parameter, 
h p 



fi = 



phri 



is the applied pressure gradient parameter, 



p. _] h 

L) = — is the inverse Darcy parameter 
k 



M 2 = — ^ — - — is the Hartmann number (Magnetic field Parameter) 
pv 



and G = ^ is the body acceleration parameter 

hn 

The shear stresses on the upper and lower walls at y = 0 and h in the non-dimensional form are calculated by 



T = jU 



+ cc. 



'3V 

dydt. 



(2.13) 



www.ijeijournal.com 



Page I 36 



MHD Flow of a Rivlin-Ericson Fluid through a Porous Medium in a Parallel Plate Channel under an 



III. Results and Discussion 

The profiles (1-6) exhibit the nature of the velocity for variations in the governing parameters S, /? , R, G, 
M and D S the Visco-elastic parameter, R the Reynolds number, /? the Applied pressure gradient parameter, M 
the Hartmann number, D~ l the inverse Darcy parameter and G the body acceleration parameter. In all these 
profiles, we observe that the velocity rapidly increases from zero value on the lower boundary y = 0 to its 
maximum value on the upper boundary y = h. For fixed R and /? the velocity enhances with increase in the 
visco-elastic parameter and this enhancement in the velocity is plotted in the figure (1) for various values of /? 
and R. For sufficiently large R the magnitude of the velocity is high comparable to smaller R and once again 
exhibits increase in trend with increase in S. The figure (2) correspond to decrease in /? for fixed R and S, we 
notice that an enhancement in the velocity for reduction in /? . Also higher values of R the magnitude of the 
velocity relatively increases for all /? and S. The figure (3) correspond to behaviour of the velocity for variation 
in R fixing /? and S, we once again notice that an increase in R enhances the velocity for fixed /? and S while its 
magnitude exhibits a slight growth for increase in the visco-elastic parameter S at all corresponding values of /? 
and R. Also the figure (4) correspond to behaviour of the velocity for variation in G fixing R, /? and S, we once 
again notice that the velocity enhances increase in for fixed values of R, ft and S. From figures (5& 6), the 

behavior of the velocity always reduces to increase the intensity of the magnetic field M, while the magnitude of 
the velocity continuously decreases with increase in the inverse Darcy parameter D" . 

The shear stresses ( t ) are evaluated for variations in R, S and ft at both the lower and upper 
boundaries. In general, the shear stress in the upper boundary is higher in magnitude compared to its value on 
the lower boundary at the corresponding sets of the parameters. From tables (1) and (2) we find on the lower and 
upper boundaries and increase in R slightly increases z for all fixed S, M, D 1 and /? Tables (3) and (4) 

represent the variation in t for increase S fixing R and /? . We notice that t enhances with S for all sets of R, 
M, D 1 and /? . Likewise, the stress z reduces with increase in /? on the both boundaries for different sets R, M, 
D" 1 and S as shown in tables (5&6). 

IV. FIGURES AND TABLES 
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Fig 1: the velocity profile for u against S with p =0.77135, R=257.8125, M=2, D 1 =1000 
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Fig 2: the velocity profile for u against p with S=0.01, R=257.8125, M=2, D 1 =1000 
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Fig 3: the velocity profile for u against R with p =0.77135, S=0.01, M=2, D 1 =1000 




Fig 4: the velocity profile for u against G with p =0.77135, R=257.8125, S=0.01, M=2, D 1 =1000 




Fig 5: the velocity profile for u against M with p =0.77135, R=257.8125, S=0.01, D 1 =1000 




Fig 6: the velocity profile for u against D 1 with p =0.77135, R=257.8125, S=0.01, M=2 
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Table-1: The shear stress at y=0 level 
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Table-2: The shear stress at y=l level 
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Table-3: The shear stress at y=0 level 
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II 
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Table-4: The shear stress at y=l level 
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Table-5: The shear stress at y=0 level 
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Table-6: The shear stress at y=l level 



V. CONCLUSION 

Based on the results of the numerical calculations, it can be concluded that: 

[1]. For sufficiently large values of R, the magnitude of the velocity is high comparable to smaller R and once 

again exhibits increase in trend with increase in S. 
[2]. To decrease in p for fixed R and S, an enhancement in the velocity for reduction in p . Also higher values 

of R the magnitude of the velocity relatively increases for all (j and S. 
[3]. The behaviour of the velocity for variation in R fixing p and S, an increase in R enhances the velocity for 

fixed f} and S while its magnitude exhibits a slight growth for increase in the visco-elastic parameter S at all 

corresponding values of p and R. 
[4]. The behaviour of the velocity for variation in G fixing R, fi and S, the velocity enhances increase in for 

fixed values of R, p and S. 

[5]. The behavior of the velocity always reduces to increase the intensity of the magnetic field M, while the 
magnitude of the velocity continuously decreases with increase in the inverse Darcy parameter D 

[6]. On the lower and upper boundaries and increase in R slightly increases r for all fixed S, M, D 1 and p . The 
variation in r for increase S fixing R, M, D" 1 and p . z enhances with S for all sets of R, M, D" 1 and p . 
The stress r reduces with increase in p on the both boundaries for different sets of R, M, D 1 and S. 
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